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Предложенный ранее [1] метод применен для вычисления всех шести 
‘независимых кинетических коэффициентов модельной дисперсной системы: 
система одинаковых, ориентированных. поперек потоков, диэлектрических 
цилиндров стонким двойным слоем. Найденные формулы учитывают поля- 
ризацию двойного слоя и являются обобщением полученных ранее формул 
‚для электроповерхностных эффектов в бесконечно разбавленной дисперсии 
на дисперсии с конечной объемной долей частиц. Установлено, что ячеечный 
метод хорошо работает. при объемной доле, ие‘превышающей 34%. 


В данной работе при помощи общего метода рассмотрения неравно- 
весных электроповерхностных явлений в концентрированных дисперси- 
ях [1] будут найдены кинетические коэффициенты (КК) следующей мо- 
дельной дисперсной системы: параллельно ориентированные бесконеч- 
ные диэлектрические цилиндры радиуса а погружены в раствор бинар-. 
ного электролита; двойной слой (ДС), окружающий частицы, предпо- 
лагается тонким. 

| ха» 1 (1) 


(х — обратный дебаевский радиус); ограничения на величину &-потен- 
циала частиц’ не накладываются. Исследуется также предел. примени+ 
мости ячеечного метода по объемной доле (на примере электроосмоти- 
ческого течения). Ссылки на формулы, работы [1] даются в виде (1), (9). 

$ 1. Конкретизация ячеечной модели. Пусть объемная доля диспер- 
сии равна %. В качестве поверхности ячейки $, выберем цилиндр радиу- 
са 6, коаксиальный с частицей, причем ‘ 


ат" (45): 


Запишем общие решения уравнений для локальных величин (1), (9) —. 
(11) для данной геометрии. Поскольку в теории поляризации тонкого, 
ДС [2] установлено, что при малых числах Пекле (Ре=Иа/О) за пре- 
делами ДС жидкость электронейтральна, а распределения потенциала 
и концентрации удовлетворяют уравнению Лапласа, мы будем искать 
потенциал в виде _ В сы 


м 


Сы т > М | 
фи = А (р) с0з 0 | (2) 
а концентрацию электролита с(г) в виде * 
с (гусь (1-Е с» (г) ) = со (1-Е с0$ 0 (@ь-+ М./р)) _ (3) 


В случае изолированной частицы М, и М. имели смысл индуцирован- 
ных дипольных моментов потенциала и концентрации, а Ёи @ — одно- 
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родных градиентов этих величин. В рамках ячеечной модели Ё, С, М. 
М.— подлежащие определению постоянные, не имеющие столь нагляд- 
ного физического смысла, как для одной частицы. 

За пределами ДС на жидкость не действуют никакие силы. Поэтому 
там скорость у, и давление 9 удовлетворяют однородным уравнениям 
Стокса и имеют вид 


ито оеучоучсь © 
и: — 9 [с, о (@)) 95, (2) о ( + е] (5) 
оне] Св 


Чтобы найти 8 постоянных, входящих в (2) —(6), необходимы 8 гра- 
ничных условий. Их выбор описан в [1], а конкретный вид таков 
1) условия прилипания (1), (12) : 


Ор в =0 (7) 
бе [р—=в= 0 (8) 


2) специфическое для конкретной гидродинамической модели усло- 
вие (1), (14) —(16) - 


ЕО 9% 1 | 
гоу =е (№ ——— =0 9 
| 7 др Г р 90 1 (9) 
в модели Кувабары, 
до о 1 99 . 
п и (= Е =) —0 10 
9 |. 1 до Г р 09/155 —_ 
в модели Хаппеля, 
Ор 0 [-ь= — 00 08 0 [рь . (11) 


в модели Мехта и Морзе; | | 
3) условия равенства нулю нормальных потоков ионов на поверх- 
ности частицы (1, 13) мы преобразуем к виду [2] 


м ть Ге ` я 
(1+ "- инету = (1— Е ) в==В)— 
а А (12) 
а 


где в = —- [ехр (= 2Ф) —  — равновесные добавки к концентрациям 
„= 


ионов в ДС (безразмерные), ф..— безразмерный потенциал в равновес- 
ном ДС (в силу условия (1) ДС будем считать локально-плоским), О*, 
2=— коэффициенты диффузии и валентности катионов и анионов (в от- 
личие от [1] в данной работе 2= обозначает абсолютную величину ва- 
лентности) ы 


= + "9 
Ги = 75 ф— а)" ар 


Полная тангенциальная скорость 9 в (12) состоит из двух частей 
бо ба 


где 0, имеет вид (5) и вызвана приложением давления, и‹,— ионоосмо- 
тическая скорость, вычисленная в $ 2 и вызванная градиентами потен- 
циала и концентрации. Каждая из скоростей у, и у, удовлетворяет крае- 
вым условиям. | 
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Обобщение (12) краевого условия теории поляризации тонкого ДС, 
полученного в [2] для изолированной частицы, возможно в силу. того, 
что ячеечная модель справедлива при 52а (см. $ 6), и при условии (1 | 
выполняется требование отсутствия перекрытия ДС соседних частиц 


(6—а)*>1 г (13) 


При выполнении условия (13) на поляризацию К данной частицы влия- 
ют только дальнодействующие поля соседних частиц. 

Пять граничных условий (7) — (12). необходимо дополнить тремя ус- 
ловиями, выражающими зафиксированные значения трех интегральных. 
величин (1), (35) —(37) через локальные поля. Чтобы получить явный 
вид этих условий и рабочие формулы для вычисления КК, необходимо: 
конкретизировать формулы (1), (24) —(26) для ячеечных потоков и (1),: 
(28) —(30) для ячеечных сил. Поскольку все указанные величины вы- 
ражаются поверхностными интегралами, для ‘вычислений достаточно: 
знать только угловую зависимость локальных величин. Поэтому форму- 

лы (14)— (19) фактически справедливы для любых поляризационных 
полей, имеющих дипольную угловую зависимость 


ЕЕ (14) 


с05 0 


| 


= И И а _ (45) 
ре р (16) 


г а Е Пов9ь (17) 
А 
= _@ 6 
и — Ь с0$ 9 | (18) 
КТе = 


В к + 19 
| Ф, $ соз 0 | Г® 


Выражение для градиента. давления (17) ` ‘можно ‘конкретизировать, ис- 
пользуя (4) — (6): к . 


— 4 5 9 й С] . С 
в [с] (|9 9 
от, .05 9 | с05 9, | Ут 0 у 
которое для ячеечных моделей Хаппеля и Мехта сводится к 
ФЕ. (20) 


95? 


Граничное условие `Кувабары (9), как отмечалось в [3], не обеспечи- 
вает «гидродинамической изоляции» ячейки, и имеется некоторая донол- 
нительная дисеипация энергии, сверх учтенной: в. ‘формуле (1), (21), за 
счет работы сил трения на`границе ячейки. Учесть этот вклад можно 
путем такого приема: создадим за пролыи ‘ячейки течение у, облада- 
ющее следующими, свойствами` 


[< при оо; Ооо ть 
И бо -ь= 9 =; Об] ь-ь==0 ) 


Теперь границу ячейки можно считать не вносящей диссипации, а толь- 
ко преобразующей внешнее течение 9 во внутреннее течение и. Добавив 
в (1), (21) и соответственно в (1), (29) член, описывающий диссипацию 

°во внешнем течении, мы.снова придем к формуле. (20) для всех ячеечных 
моделей (с точностью до членов 9(=)). 


554 о а ть а р 


© 


Величина | Модель Кувабары Модель Хаппеля Модель Мехта 


Коэффициенты поля 
скорости 99 


7 —21, —21, а 
, 2 
9 аз —2 т 42 | 2-я |а) пу 4—2) 
не |-% | Е р 
в= 1 яп 0 бр | ад —\) 799 АР. —в0 #—”) 
Поправочный множи- 4 1-8. - 1 
тель Н (<) в (1), (35) У. АНЯ 
Коэффициент поля ско- 
рости 91 й . ? ; , 
— ме с: 
С. Ио (41—72) | | ея ) Уз (4-- =) 


Формулы (14) —(16) и (18) —(20) завершают `конкретизацию ячееч- 
ной модели. 

Выпишем для справок наиболее важные гидродинамические величи- 
ны для рассматриваемых моделей (таблица): 

$ 2. Учет ионоосмотических поправок в течении. Вычислим ионоос- 
мотическое (т. е. электро- и капилярноосмотическое) течение у, которое 
необходимо для конкретизации уравнений (12). Оно имеет вид 


У==У, НУ» 95==0 


р 4 №, о ге. ье 9 (21 


где рз = — ЮТ (2*с* — ат ф, и = ЮТсх, ей —2 у) Фы— про- 


екции объемной силы у.— Решение вида (4), (5), Е кото- 
рого определяются из условий (7) —(11) на скорость 9; и Оно [в-=ь=2 0. 

Коэффициент С., входящий в У; и определяющий. поправку к гидро- 
динамической проницаемости, выписан в таблице: (У = — 0 | р=о/31п 6). 
Формула (21) по существу совпадает с аналогичной формулой для ионо- 
осмотического течения, выведенной в [4]. Можно пора Зоват» (21) к 
виду 


ВТе 


|(-2 +) тнг го) + 


зв = 


+(5+ ат г |. 9 (22) 


где Г! (0) = й уе (р — х)ах. Как видно из (22), бы экспоненциально убы- 


вает за пределами ДС, а на поверхности частёцы принимает значение — 
У, зш 0. 9. медленно убывает от значения ТУ, $10 на поверхности час- 
тицы до своего значения на поверхности ячейки. При вычислении вклада 
скорости оне в интегральный конвективный поток в (12) достаточно огра- 
ничиться старшими по обратному ха членами, поскольку они имеют по-_ 
рядок Ке{. Поэтому интеграл от 0. надо брать точно, а в интеграле от 
о.. скорость можно разложить в ряд Тейлора на`поверхности частицы 
и ограничиться первым ненулевым членом, Т. е. 9». (а), который не зави- 
сит от гидродинамической ячвечной модели. Таким образом, в рассмат- 
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риваемом приближении (тонкий ДС, учет поляризации с точностью до 
О (1/ха)) все ячеечные модели дают одно и то. же значение ионоосмо- 
тической поправки. При вычислении интеграла от 9, также достаточно 


до 


| [] 
ограничиться первым неисчезающим членом разложения, т. е. - о 


р: 
С учетом сказанного (12) можно переписать в виде 


М, (1 -- В*) — а?6 (1 — В=) 4 2==Мо(1 + А*) + 242 (1— А*) = 


се 
2= 


а Е _ 
= ЕВА (23) 


где 


ВТ К=. В=— т + ВТег= КЕ 
РЫСИ бк КЕ 


т2=а а тр=а 


К= = — [1% (х) ах { (х— 1 (27—20 0)@ 


ке фа [ео Оо-тОи 


х 


$ 3. Электропроводность дисперсии и диффузионный потенциал. За- 
дадимся значениями «ячеечных» градиентов } 


Ф,=0 (24) 
Ф,=сопз{` | | (25) 
Ф,=0^ > (26) 


тогда в соответствии с уравнениями ‘переноса (1), (31), электропровод- 
ность 9, равна 


ги 
| “=, 
а коэффициент .:, характеризующий диффузионный потенциал, равен 
дить, 
| аа | 
Уравнения (25), (26) в соответствии с (18), (19) и (2), (3) приводят к 
Е М, 
а = — г (27) 


Учет условия (24) в уравнениях (23) вместо Г.=0 приводит к поправ- 
кам, в ха раз меньшим учтенных. Поэтому мы положим Г,=0 и решим 
систему уравнений (23), (27), (28). Подставив поляризационные поля 
последовательно в (1), (8); (1), (22), (23); (15), (16), получим 


Фа р РТВ) +" ее] 9 


\1-- А 1, 
— Во. — ээ [1-9 2+ 27) 4рк, 
* — кт [= м в | р 


где } 

ГВ, = [2*(1—А+) /М++- = (1 А-)/М-М=*(1+А+) /М++=-(1+А-)/М-] 
Ре. Г, Еы-Р.Р- 2 к+к-|/. | 

ь-| ее (>=) р АНИ АТЧНА 


а а? 21а / Р+р- 
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аа 
при с.-—>со; М=—1- 12; Ю,->0; Ю—1 
и й 
Е но 
9 (гр Нар) = (“р+*- г0-) (1—2) (31) 


Как видно из (31), электропроводность (28) удовлетворяет требованию 
(1), (32). Это означает, что интегральные потоки /[, и силы Ф, совпада- 
ют с «ячеечными» (для #=2,3). Таким образом, равенства (1), (36), 
(37) доказаны. 

$ 4. Гидродинамические эффекты. Определение поправочного мно- 
жителя Н (). Зададимся значениями градиентов 

Ф, =сопзЕ; Ф,=0; Ф,=0 

Тогда можно вычислить следующие КК: гидродинамическую проницае- 
мость 9, =— Г,/Ф,; ток течения (электроосмос) ,„=—Г,/Ф;: капилляр- 
ный осмос 9, =—Г,/Ф.. Действуя аналогично 6 3, получим 


Зичи О-ВА — Фет! ТР ВНЕ Во] " 
и а ‚ 


4 
(32) 
где 
м 
ми; В; ГЕ) 
в = 2% 980 [++ 07) В, + (2*—Р^ Ви (33) 
21 1 (9) 
— -%_ 98% [(р+ _р- а 34 
2 Г |" В+ (++) (94) 
где 
СГ Г. '|/| 2М Му 
ие р+м+  р-м- И м+ и 


` Г+ г: - 
Ва а = 
р+м" РМ: РМ 2М; 
М.== 1 + (1—4?) А= 


При с.—с°, ха—0, Ке1->0 условия (1), (3), (4) выполняются, при этом 
М:=—1-5?, 


рев] 


— рва я 
За _ = ов (ра (35) 
ба 4ллК 21 (9). (1 %) 


Структурно-зависимый множитель в (35) не равен единице. Однако если 
мы введем поправочный множитель Н (о) в соответствии с (1), (35) и 
приравняем его | 
а _ 98 (у) : | 
Н (@) = — (36) 
©) 2 1) а—\") 


то формула Смолуховского будет выполняться для любой ячеечной мо- 
дели. Значения множителя Н (а) приведены в таблице. Таким образом, 
удовлетворен критерий (1), (33). В соответствии (1), (31); (1), (34); 
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(1), (35) интегральные КК Г.» связаны с ячеечными КК т формула- 


ми (37 | 
Га а/Н? (а); Гь= Зы Н (а) (#=1--3) 

[ит = „т (п, т=2,3) (37) 
| $ 5. Диффузионная проницаемость. Зададимся значениями гради- 
ентов ы № А | . 

Ф,=0, Ф.=0,. Ф,=соп$ 
Тогда диффузионная проницаемость .. равна 


: 5? а —7 з/ Ф, 
Действуя аналогично $ 3, получим | | 
__1 0+ р 1—9 \ _ 2+— р 22 38 
® ЭТ | ЕЕ К 1 в, а, 05 


где 


тв 1-1 ака с ав 
ыы [ АМ 2: |, [ Ам ЕМ 


1—8* 1 8-|/| 2 М: ги: | 
В — 1+8+ 1В- я Е 1+ 8- 
Таким образом, формулы (29), (30), (32) —(34), (38) дают значения 


«ячеечных» КК, а формулы (36), (37) устанавливают их связь с интег- 
ральными КК. 


5,% 


Зависимость  погрешно- 

сти ячеечной модели при 

описании  электроосмоса 

от объемной доли частиц 
- в дисперсии 


ШИ?” 90. 90 70 =, % 
! 

$ 6. О пределе применимости ячеечного метода. Интересная возмож- 
ность оценить точность ячеечного метода возникает при сравнении поля 
скорости электроосмоса, найденного. на основе ячеечной модели, и не- 
зависимого приближенного выражения, построенного с помощью кон- 
формного отображения [5]. Последний метод является регулярным и в 
принципе позволяет получать результат с любой точностью. 

В указанной книге приведена таблица значений 9. |,-. для обтека- 
ния решетки цилиндров при различных значениях 1; чтобы получить. 
аналогичную величину в ячеечном методе, необходимо решить систему 
уравнений (7), (9)—(11), (24) и | 

| бр |р=ь=1: С0$ 0 ` 

Легко видеть, что ) 
21; $т 9. 
1—9. 


9 —а = — (39} 
Возьмем в качестве меры погрешности р 
п/з. ма а : 
6() = | (19 [-: =. А 1 5).40, [ [95| „99 | 
Г) | т. ое а 
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где 9е|р=« берется по (39), а о |›.— из табл. 11 в работе [5]. График 
5(«) приведен на рисунке. Как видно, погрешность ячеечной модели мо- 
нотонно возрастает с ростом объемной доли. Если считать погрешность 
в 10% допустимой, то можно сказать, что ячеечная модель хорошо опи: 
сывает поле скорости электроосмоса вплоть до «= 0.34. 
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Зишшагу 


Тре шеёо4 ргорозеё еагНег [1] Ваз Веки: аррНей 4 сайсШа{е аП эх ш4ерепйет 
стейс сое сет оЁ а то4е| 41зрегзе зуз4ет:` {Ве зузёет “о! ‘еана; {тапзуегзе!у оцещей 
Угеатз, Феесийс суНп4егз %УЙВ а Нил доиЫе 1ауег. Тне омти!аз дейуе4 4аКе {о ас- 
сои ро!апхаНоп оЁ {пе донЫе 1ауег ап@ аге а сепегаН2аНоп о{ 4Не югти!аг оМатей 
еагНег {юг е еесёгози{асе еНесё$ ш ап ш—пйеу 4Ищей @1зрегзюп оп Не 41зрегз!оп 
Вауше Ше НаНе БшК гасНоп оЁ рагНс!ез. И Ваз Бееп езфаБИзНе@ 4па{ {Пе сеПёаг те{- 
Воф ре{огилз уеИ ИН а БШ {гасНоп {На{ 4оез по{ ехсееё 34%. 
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